Multiple-charge transfer and trapping in DNA dimers 
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We investigate the charge transfer characteristics of one and two excess charges in a DNA base-pair 
dimer using a model Hamiltonian approach. The electron part comprises diagonal and off-diagonal 
Coulomb matrix elements such a correlated hopping and the bond-bond interaction, which were 
recently calculated by Starikov [E. B. Starikov, Phil. Mag. Lett. 83, 699 (2003)] for different 
DNA dimers. The electronic degrees of freedom are coupled to an ohmic or a super-ohmic bath 
serving as dissipative environment. We employ the numerical renormalization group method in 
the nuclear tunneling regime and compare the results to Marcus theory for the thermal activation 
regime. For realistic parameters, the rate that at least one charge is transferred from the donor 
to the acceptor in the subspace of two excess electrons significantly exceeds the rate in the single 
charge sector. Moreover, the dynamics is strongly influenced by the Coulomb matrix elements. We 
find sequential and pair transfer as well as a regime where both charges remain self-trapped. The 
transfer rate reaches its maximum when the difference of the on-site and inter-site Coulomb matrix 
element is equal to the reorganization energy which is the case in a GC-GC dimer. Charge transfer is 
completely suppressed for two excess electrons in AT- AT in an ohmic bath and replaced by damped 
coherent electron-pair oscillations in a super-ohmic bath. A finite bond-bond interaction W alters 
the transfer rate: it increases as function of W when the effective Coulomb repulsion exceeds the 
reorganization energy (inverted regime) and decreases for smaller Coulomb repulsion. 

PACS numbers: 71.27.a, 34.70.e, 82.39.Jn 
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I. INTRODUCTION 

Understanding charge transfer dynamics in DNA is 
of fundamental importance for biochemistry and molec- 
ular electronics It can generate chemical reactions 
over long distances and has biological relevance in the 
formation of oxidative damage. Charge transfer in 
DNA also plays an important role in anti-cancer drugs, 6 
electrochemical readout of micro arrays and molecular- 
electronic devices. Whereas DNA conductivity is still dis- 
puted, there is no doubt that electron transfer through 
DNA double strands takes placed It is generally be- 
lieved that the electron or hole migrates by super- 
exchange tunneling or multistep hopping in a donor- 
bridge-acceptor system depending on the energy differ- 
ence between the bridge and donor/acceptor— and the 
fluctuation of the environment 

Most of the current understanding of electron trans- 
fer is derived from a single-excess charge where only two 
states are needed to model a donor-acceptor system. At 
high temperatures the transfer is governed by thermal ac- 
tivation. The rate can be described in the non-adiabatic 
limit by the semi-classical Marcus theory^ and is deter- 
mined by the hopping integral A, the energy difference 
between product and reactant states e, the reorganization 
energy E a and the temperature T. 

At low temperatures, however, the transfer rate does 
not vanish but rather remains finite and temperature in- 
dependent: the transfer is governed by quantum fluctua- 
tions and involves nuclear tunneling^ In this regime the 
transfer is modeled by an electronic two-level systerr&i£ 



coupled to a dissipative environment as in the spin-boson 
model 14 which has to be treated quantum mechanically. 
One can identify three regions of dynamics depending 
on the coupling to the bosonic bath: coherent (damped 
oscillations) , incoherent (exponential relaxation) transfer 
and self-trapping. 

However, in many cases charge transfer involves multi- 
ple excess charges which are correlated i 1 ^— Their trans- 
fer characteristics is changed drastically in comparison to 
the case of a single excess charge and requires a descrip- 
tion by more than two states.— ~— Two charges can mi- 
grate as pairs, sequentially or can be self-trapped. Their 
dynamics are governed not only by the coupling to a dis- 
sipative bath but also by the local and non-local Coulomb 
interactions (U and V, respectively. ) 16 i 17 i 21 

In this paper we study the charge transfer dynamics 
of two electrons or two holes in different DNA base- 
pair dimers based on Guanine/Cytosine (GC) or Ade- 
nine/Thymine (AT) pairs. We have used the base-pair 
specific Coulomb matrix elements as calculated with ab- 
initio methods^ We assume that two charges are ini- 
tially prepared on one base-pair monomer (donor) - this 
could be achieved, e.g., by voltage gating. Then, the 
transfer of those two charges occurs either in a concerted 
manner, where an electron pair or hole pair is transferred 
between the donor (D) and the acceptor (A), or both 
charges migrate sequentially but correlated. This trans- 
fer of the first electron facilitates or hinders the transfer 
of the second one due to the electronic correlations. We 
find that depending on the Coulomb matrix elements a 
second charge can lead even to self-trapping or, in the 
other extreme, to partially activationless transfer. 
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Recently, conductance measurements of single DNA 
molecules were performed in aqueous solution^ The so- 
lution preserves the native confirmation and the trans- 
port properties can be associated with charge trans- 
fer reactions at least for weak coupling to the leads i24 
At voltages equal to the on-site Coulomb interaction 
we expect that more than one excess charge is rele- 
vant for the charge transport through a molecular de- 
vice. Consequently, the transport characteristics is al- 
tered considerably^ Furthermore, two-electron transfer 
in Guanine due to strong oxidants has been discussed in 
connections to in vivo lesions^ 

Our detailed study of the base-pair specific transfer 
characteristics of multiple excess charges for realistic sys- 
tems reveals the connection between the microscopic pa- 
rameters and the transfer efficiency. Previously, only an 
electronic model with density-density interactions cou- 
pled to an ohmic bath has been investigated^ In this 
paper we extend the calculations to a more general elec- 
tronic model which includes an additional bond-bond in- 
teraction (pair-transfer) W previously neglected as well 
as a bond-charge interaction X (correlated hopping.) 
These matrix elements have been calculated for differ- 
ent DNA base-pairs^ 2 - using ab-initio methods and serve 
as input for our calculations. 

We have found that the additional bond-bond interac- 
tion W favors concerted electron transfer from the donor 
to the acceptor but can also lead to a change of the 
ground state of the dimer. The parameter W has a pro- 
found impact on the ground state of the electronic sub- 
system with two excess charges: it lifts the degeneracy 
of singlet and the triplet state. The latter becomes the 
ground state for W > and a positive effective Coulomb 
matrix element. Since the Hamiltonian considered here 
is spin-conserving, the system can only decay into an 
metastable steady state in this regime starting with two 
excess charges on the donor in a spin singlet state. 

Another important application of charge transfer are 
DNA biosensors. During the sensing process, electron 
transfer is used to produce electronic signals^ In gen- 
eral the signal magnitude depends on the efficiency of the 
electron transfer through the DNA as well as the number 
of electrons. It has been suggested that the transfer of 
multiple electrons would enhance those signals i21 How- 
ever, our work shows that the transfer rate for a single 
electron in systems with two excess charges depends on 
the environment as well as the Coulomb matrix elements: 
it may be much larger or much smaller than the transfer 
rate of a single excess charge depending on the micro- 
scopic parameters. 

For temperatures lower than the activation energy 
the bath must be treated quantum-mechanically to cal- 
culate the time-dependent population probabilities and 
the transfer rates. In this regime, we employ the 
time-dependent numerical renormalization group (TD- 
NRG) approac h 17 ! 28 " — whereas in the thermal activation 
regime we compute the transfer rates using the Marcus 
theory 



The paper is organized as follows. In section III Al we 
introduce the microscopic model as well as both meth- 
ods for high and low temperature regimes. We present a 
detailed description of the single and multiple charge dy- 
namics in Sec. IIIII and Sec. IIV1 respectively. A summary 
and conclusions of our results is given in Sec. [VJ 



II. MODEL AND METHODS 

A. The model 

We consider a base-pair dimer represented by a two- 
site system as donor (D) and acceptor (A) with a LUMO 
(lowest unoccupied molecular orbital) and a HOMO 
(highest occupied molecular orbital) coupled to an aque- 
ous environment (bosonic bath) as sketched in Fig. [T] 
The Hamiltonian of the coupled system consists of three 
parts 



H — H e i + Hb + H ( 



l-b 



(1) 



where H e i denotes the electronic part of the Hamilto- 
nian, Hb the bath and H e i— b the coupling between both 
subsystems. Each base pair monomer is modelled as a 
single site with one orbital (either LUMO or HOMO de- 
pending on excess holes or electrons). The energy gap 
between HOMO and LUMO is about AeV £ H e i has been 
derived in Ref . [22| and corresponds to a general two-site 
extended Hubbard model: 

H el = + J2 A '^W ( 2 ) 



i,j,k,l,a,<r' 



where and c ia denote annihilation and creation opera- 
tors for fermions with spin a in an orbital localized at site 
i. The first and second term are the on-site energy and 
the transfer integral or single-charge hopping for near- 
est neighbors, respectively. The latter defines the kinetic 

kl > account for the 



energy. The matrix elements {ij 

Coulomb repulsion between two electrons. The dominant 
interactions in the DNA dimcrs — are 



Ui 



(3) 



The first two matrix elements denote the on-site and 
nearest neighbor electron-electron interaction, respec- 
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FIG. 1: (Color online) Schematic view of the DNA dimer. 
Each base pair is approximated by a single site with a HOMO 
and LUMO. The electronic degrees of freedom are coupled to 
the environment comprising internal vibrations and the sol- 
vent dynamics. The transfer characteristics of multiple excess 
charges strongly depends on the base-pair specific Coulomb 
interactions. 



tively. The parameter X labels the bond-charge inter- 
action and is often called correlated hopping since it de- 
pends on the density of the nearest neighbor. The last 
term is the bond-bond interaction which induces the pair 
hopping and the exchange interaction. 

For the donor-acceptor pair, we arrive at the following 
Hamiltonian: 



H e l = ^ e i n ia- - A ^ [ C Da C Aa + C Aa C D<y 
a;i=D,A a 

+ U n ia ni-a + V n D<? n A<j' 

i=D/A,a o",er' 

+ X y (cpgCAcr + C^Cpa) {n D - a + n A - a ) 
a 

+ W Yl ( C D* c A-a C D-*CA* + c\ a C X D _ a C A - a C Da 

(7 

+ £Da C D-a C A-*C A * + c\ a c\_ a C D - a C Da ^ (4) 

with the number operator = ct c- a . The last two 
terms correspond to hopping of electron pairs. In our 
previous studyjii only the matrix elements U and V have 
been taken into account. Here, we will investigate (i) the 
influence of the additional matrix elements X and W 
onto the transfer dynamics using (ii) realistic values for 
U, V, X, W as well as for the coupling to the environment 
as will be discussed in the following. 

The bosonic bath comprises the internal vibrations and 
solvent dynamics and leads to dephasing and dissipation 
of the electronic system. The Hamiltonian Hb 



(5) 
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FIG. 2: Schematic view of the relevant states: The ionic states 
j t-U 0) and 10,4/f) as well as the covalent singlet state 
They are connected by the renormalized hopping A — X as 
well as pair hopping W. |^t) is not connected to the other 
states. The energy difference between the ionic states and 
\^ s ) is UcS — W. Only the ionic states are coupled to the 
bosonic bath. 
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and give rise to an attractive contribution to the on-site 
Coulomb interactions. The renormalized interaction is 
defined by 



U cS = U - V - 



A 



(7) 



which determines the energy difference between the ionic 
and covalent states for W = and includes the influence 
of the environment. 

The Hamiltonian of the isolated dimer conserves the 
total charge Q = 4 — n, where n is the total number of 
electrons. We consider the case Q = 2 (two holes (2h) in 
the HOMO), Q = -2 (two electrons (2e) in the LUMO) 
as well as Q — 1 and Q — — 1 (single hole in the HOMO 
and single electron in the LUMO, respectively). 

Furthermore, the square of the total spin as well as 
its z-component S z are conserved by the Hamiltonian 
Eq. ©. The subspace of two holes/electrons (|Q| = 2) 
and S z — (for \S Z \ = 1 no transfer is possible) is 
spanned by 4 states, two ionic states with two elec- 
trons/holes on the donor | fl, 0) or on the acceptor |0, t>L) 
and the two covalent states 



l* s > = ^(It4)-U,t» 
l*t) - ^(|t,+) + U,t» 



(8) 



models the free bosonic bath. The displacements of 
the polaronic bath couple to the change of the electron 



The latter states correspond to the singlet and triplet 
states, respectively. \^t) is an eigenstate of H e \ with 
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energy E t — —W and is decoupled from the other three 
st ates displayed in F ig. [2] It is the ground state for 4W > 
y/U 2 + 16(A-X) 2 -U (a = 0). 

In the subspace of \Q\ — 2, S z — 0, the total number 
of electrons per spin, nu a + nAm is constant, and the 
correlated hopping amplitude X renormalizes the single 
particle hopping to A — » A — X. 

For one electron or hole (\Q\ = 1) and \S Z \ = 1/2 
the model is equivalent to the spin-boson model. The 
charge is able to migrate between the two states |», 0) and 
|0, •) with one electron/hole on the donor and acceptor, 
respectively. 

The coupling between electronic subsystem and 
bosonic bath is specified by the bath spectral 
function! 4 -^ 



(9) 



J(uS) characterizes the bath and the system-bath cou- 
pling and reflects intramolecular vibrations and the sol- 
vent (water). The spectral function can be calculated 
by molecular dynamics simulations^ for a two-sphere 
DA model embedded in an aqueous medium and shows 
a power law behavior at low frequencies. The asymp- 
totic low-temperature behavior is determined by the low 
energy part of the spectrum. Discarding high-energy de- 
tails, we use the standard parametrizatio n 14 ' 36 



J(u) 







< UJ < LJ C , 

otherwise , 



(10) 



where a is a dimensionless constant which characterizes 
the dissipation strength and u c is the cutoff frequency of 
the spectral function. 

We investigate an ohmic bath (s = 1) and s — 2 as an 
example for the super-ohmic dissipation^ As an realistic 
estimate for the cut-off frequency oj c we use 0.5eV which 
correspond to the frequency of the OH stretch vibrations 
at about SmQcmr 1 ^ 

The classical reorganization energy E a for the single- 
electron transfer |»,0) — > |0, •), which measures the en- 
ergy relaxation following a sudden electronic transition, 



„ w « n Jo 



did 



J(uj) 2auj c 



(11) 



has been determined 3 for DNA in water as approximately 
E a = 0.5eV. Therefore, we have used a = 0.5 and a = 1 
for the ohmic bath and super-ohmic (s = 2) bath, re- 
spectively, for modeling the coupling between the DNA 
dimers and the environment. The reorganization energy 
for the transfer | fi, 0) — > \9 S ) in the sector of two ex- 
cess electrons is given by E a \ = 4E a = 2eV and for the 
transfer |tl,0) -> |0,t4) given by E a2 = 4E Ql = 8eV. 

The single-charge hopping A as well as pair-hopping 
W are renormalized by the bosonic bath. Because only 
the ionic states are coupled to the bath, the pair-hopping 
is much stronger reduced than A which connects ionic 



10 



LU 




-4 -2 





2 


a 








X 








8 


W 1 












J T40^ 










LU 










w 





1 — i 


1 


1 — 








4 -2 




X 


2 


t 


. 



FIG. 3: Schematic potential surfaces E(x) for the lowest 
states of the model in the one-electron (upper panel) and 
two-electron subspace (lower panel) for the displacement x 
in arbitrary units. The energy difference between the ionic 
states (| t4-j0), |0, t4-)) an d covalent singlet state corre- 
sponds to the difference between the renormalized Coulomb 
interaction U c s and bond-bond interaction W. 



and covalent states. The exchange interaction, however, 
which is also proportional to W remains unaltered by 
the bath since it does not change the occupancy of the 
orbitals. 

In the following sections, the four time-dependent pop- 
ulation probabilities, defined as the expectation values of 
the operators 

n D = h,0)(0,.| 

d D = in,o)(o,m 
d A = io,nxn,oi 

tida = |t,4>U,t| + U,tXt,l|, (12) 

are calculated. They measure the single and double oc- 
cupancy on the donor and the acceptor site (ionic states) 
as well as the equally populated states (covalent states). 
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B. The Marcus theory 

The Marcus theory of charge transfer is connected to 
the quantum mechanical treatment of the bosonic bath 
in Hamiltonian Jl) by a mean-field decoupling of the 
bosonic bath modes from the electronic degrees of free- 
dom. The local displacement x = 2^ n A n (6t + b n ) is 
replaced by its expectation value x, which is also called 
the reaction coordinate. The resulting potential surfaces 
E(x) for the two states D~ A = |»,0) and DA~ = |0,») 
in the |Q| = 1 sector are plotted schematically in the 
upper panel of Fig. [3] The potential surfaces E(x) of 
the four states in the \Q\ = 2 sector, D 2 ~ A = | 0), 
DA 2 - = |0,tJ.) and D~A~~ = |* a ) or |* t ) are depicted 
in the lower panel of Fig. [3] 

The single-electron transfer from D~ A to DA~ is de- 
termined by the hopping A and the coupling constant 
a (or the reorganization energy E a ) . For \Q\ = 2, the 
transfer additionally depends on A — X and W. More- 
over, the effective bath coupling constant is four times 
larger due to a factor 2 in the coupling term. As a conse- 
quence one would expect a slower transfer compared to 
the \Q\ = 1 case. However, the transfer of two charges 
is strongly affected by the energy difference between the 
ionic and covalent singlet states U e g — W. The transfer 
of one charge can be even activationless if the difference 
U c s — W is equal to the reorganization energy E a i as 
indicated in the lower panel of Fig. [3] at finite values of 
x. 

To address the thermal activation regime where the 
temperatures are larger than the activation energy (T > 

^ Eal 4b" W ^ )■> the dissipative environment is treated 
classically, and the rates are obtained from Fermi's golden 
rule. Generally the Marcus rate is given by 

h^f = — \(f\H DA \i)\ 2 F(AE if ), (13) 
where in the Franck- Condon integral 

F(AEif) = =e *Eirf kB T ( 14 ) 

\j^ET S k B T 

the energy difference between state i and /, AEif, and 
the corresponding reorganization energy E l £*f - as de- 
scribed after Eq. (jTTj) - enters. The population probabil- 
ities are determined by kinetic equations 

P s = S ' k r ^ s P r — S ' fc s _j. r P s , (15) 

r r 

where P s is the probability that the system is found in 
state s and fc r ->s is the Marcus rate of the transition from 
state r to state s. 



C. Time-dependent numerical renormalization 
group approach 

At low temperatures (T « ^^') the bosonic 
bath must be treated quantum mechanically In this 
regime, coherent quantum oscillations may be present, 
and we have calculated the time-dependent population 
probabilities using the TD-NRG JI^^ 

For the models studied in this work, the impurity cou- 
ples to a bosonic bath, for which the standard NRG ap- 
proach has to be modified as described in Refs. I3h - l33l 
The basic features of the bosonic NRG are as follows: (i) 
the logarithmic discretization of the bath spectral func- 
tion J(ui) in intervals [A - ( ,l+1 ) , A - ™], with n = 0, 1, oo 
and A > 1 the NRG discretization parameter (all the re- 
sults shown in this paper are calculated using A = 2); 
within each of these intervals only one bosonic degree of 
freedom is retained as a representative of the continuous 
set of degrees of freedom, (ii) The mapping of the result- 
ing Hamiltonian onto a semi-infinite chain, (iii) The it- 
erative diagonalization of the chain-Hamiltonian via suc- 
cessively adding one site to the chain. 

For the results shown here, we use iV& = 8 bosonic ba- 
sis states for each chain link and retain the N s = 100 
lowest eigenstates after each iteration. The set of all dis- 
carded states eliminated during the NRG procedure form 
a complete basis of the semi-infinite chain, which is also 
an approximate eigenbasis of the Hamiltonian. In one 
NRG run we compute the initial density matrix , 17 ' 28 ! 29 in 
the second run the approximate eigenbasis of the Hamil- 
tonian Eq. (Q| . Since the Hamiltonian usually does not 
commute with the initial density matrix, the time evolu- 
tion of the density matrix governs the time dependence 
of all expectation values. 

The real-time dynamics of the population probabilities 
has been calculated by summing over all contributions 
from the sequence of reduced density matrices which con- 
tain the information of decoherence and dissipation at a 
given energy scale w c A~™. To simulate a bath continuum 
we have averaged over different bath discretizations (we 
used eight different realizations) . For more details on the 
approach see Refs. [281129} . 



III. SUBSPACE OF A SINGLE CHARGE 

Before analyzing the transfer of multiple excess charges 
we briefly review the single-charge transfer. The two- 
state system is prepared such that the donor is occupied 
by one electron/hole at t = 0. The transfer of a single 
charge is not effected by the Coulomb matrix elements 
and Eq. ([3]) reduces to the first two terms which is equiv- 
alent to the spin-boson model. We use a typical hopping 
matrix elementA = 0.4cj c = Q.2eV as obtained by ab- 
inito methodsiS 2 . and a = 0.5 as estimated above. 

The population of the donor no (t) calculated with the 
TD-NRG for T ^ is shown in Fig. H After about 200 
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FIG. 4: (color online) Time-dependent population probability 
of the donor site, nn(t) in the subspace of a single charge. The 
system is prepared such that at time t = the charge occupies 
the donor. The results are shown for E a = 0.5eV (blue), a 
decoupled charge (dotted line) for a = 0. The inset shows a 
comparison with simple second-order perturbation theory in 
a which is only valid for very short time scales. 

fs the system relaxes to its ground state Ud = h- The 
frequency of the oscillations increased due to the strong 
coupling to the bosonic bath, which is counterintuitive 
at first. 

Usually, a reduction of the frequency in a damped har- 
monic oscillator is expected due to friction. This has 
also been found in the spin-boson model in the limits of 
A/uj c — > 0, when the bare oscillation frequency is much 
smaller than the bath cutoff energy. In that limit, all 
coherent oscillations are suppressed once a — > 0.5. 

By comparing the result to the dynamics of a decou- 
pled spin (a — 0, dotted line) we find indeed a slower 
spin-decay at ultra-short time scales of t < 4fs - see inset 
of Fig. |U However, the realistic DNA parameters yield 
2 A ss lu c : the bare frequency of oscillation is of the or- 
der of the cutoff of the bosonic bath. Now, band edge 
effects become important, and the strong coupling to the 
bath (a — 0.5) lead to a level splitting in the first RG 
steps (level repulsion). Some local excitation energies 
are shifted above the cutoff lu c ; other contributions are 
renormalized further down to low frequencies in each RG 
step. This is reflected in the Fourier transformation of 
P{t) - not shown here - which peaks at a frequency of 
approximatly 1.1 3w c ; however, it also comprises a second 
broad contribution which has its maximum at a reduced 
frequency of ui ~ 0.18cj c . 

The analytic solution obtained in second-order pertur- 
bation theory in the bath-coupling \ q agrees very well 
with the ultra-short dynamics of TD-NRG. This is de- 
picted in the inset of Fig. 2) the perturbative and the 
TD-NRG solution are identical for t < 3fs for ohmic and 
super-ohmic baths - not shown here. In fact the pertur- 
bative solution shows even a stronger renormalization of 
the frequencies than the TD-NRG. However, the solution 



Dimer 


UeS 


AT-AT (2h) 


0.1 eV 


AT-AT (2e) 


-1.0 eV 


GC-GC (2h) 


2.0 eV 


GC-GC (2e) 


1.0 eV 



TABLE I: Coulomb interaction parameter taken from 
Ref. [H]. The other parameters are A = 0.2eV, X = O.leV, 
L0 C = 0.5eV, E a = 0.5, E al = 2eV and E al = 8eV. We vary 
W between -0.25eV and 0.25eV. 



becomes less reliable for t > 6fs due to the strong cou- 
pling effects ignored in the perturbation theory, which is 
only valid at very short time scales. The origin of this 
enhancement can be traced to a factor A/(cj 2 — A 2 ) in 
the analytic integral kernel which causes a sign change 
of the contribution for cj < A to cj > A. The position 
of A relative to the cutoff w c and the power s of J(ui) 
determines the sign of the frequency renormalization. 

The spin-boson model is usually only investigated in 
the limit A/w c — > 0. The DNA dimer system, however, 
appears to be in the complementary regime where 2A — > 
uj c . In the NRG as well as in the analytical solution, we 
find coherent oscillations which prevail well above a = 1 
(not shown here). 



IV. SUBSPACE OF TWO CHARGES 

To investigate the influence of the diagonal and the ad- 
ditional off-diagonal Coulomb interactions on the trans- 
fer of two charges, we initially prepare the system such 
that two electrons/holes are localized on the donor at 
time t — 0. This can be achieved experimentally, e.g., 
through radiolytic or photolytic excitation, strong oxida- 
tion or voltage gating. 

The hopping matrix element A is renormalized to A — 
X in the presence of a finite correlated hopping matrix 
element X in the subspace of two charges. The base- 
pair specific Coulomb interaction used in the following 
are stated in table I and are taken from Ref. [22j]. Since 
we are mainly interested in the influence of the Coulomb 
interactions onto the transfer rates, we set the hopping 
matrix element to a typical valued A = 0.2eV and a 
correlated hopping parameter X = O.leV for all base- 
pairs^ 2 - the dynamics of a single charge is governed by 
the energy scale of O.leV. 

It has been argued 2 ^ that the bond-bond interaction 
W is either small or negligible. Therefore, we have in- 
vestigated the influence of W onto the electron trans- 
fer dynamics by varying W between W = —0.25eV and 
W = 0.25eV. In the next section, we begin with W = 0. 
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FIG. 5: (Color online) Rate koA as defined in Eq. (|16[) at high 
temperatures (black line) and low temperatures obtained from 
NRG (red points) , as a function of U e g . The region of sequen- 
tial and concerted pair transfer is indicated. For U e g > 6eV 
the pair transfer is followed by a very slow transfer to the co- 
valent state. At low temperatures three regions are present: 
Self trapping for an ohmic bath (left of the dotted line), se- 
quential transfer (white region) and pair transfer combined 
with slow sequential transfer (shaded region for t7 e ff > 0). 
The low-temperature rates were obtained by fitting the pop- 
ulation probabilities to an exponential function. The points 
represent the rates for (1) AT-AT(2h), (2) GC-GC (2e) and 
(3) GC-GC (2h). 



A. Zero bond-bond interaction 

The Marcus theory is applicable in the temperature 
regime where the bosonic bath can be replaced by a sin- 
gle reaction coordinate, and the physics is governed by 
thermal activation. The theory predicts an exponential 
dependency of the transfer rates on the energy difference 
U e f[ — E a \. The concerted transfer of both electrons be- 
comes relevant once the transfer rate between the ionic 
states exceeds the rate between ionic and covalent states. 
The latter regime, however, is excluded by the parame- 
ters stated in tab. [Jin the DNA dimers considered in this 
paper. 

The rate hp a that at least one electron is transferred 
from the donor to the acceptor in the thermal activation 
regime is defined as 



^DA 



(16) 



and is plotted versus U e g as a solid line in Fig. [5] The 
maximum at £/ ff = E a \ marks the optimal or activation- 
less regime and separates the normal regime, f7 c ff < E a \, 
from the inverted regime, U c g > E a \. For U e g < —2eV 
and U c s > 6eV the pair transfer is the main process. For 
U Q g > it is accompanied by a slow sequential transfer 
to the covalent state. 

Using the parameters of tab. HI we note that the trans- 



FIG. 6: (Color online). Low-temperature population proba- 
bility of the doubly occupied donor P(t) = do(t) (starts at 
one) and the covalent states P(t) = nr>A(t) (starts at zero) for 
AT- AT (2h) (red dotted line), GC-GC (2e) (black solid line) 
and for GC-GC (2h) (blue dashed line) as a function of time 
for an ohmic bath (s = 1) and W = using the TD-NRG. 
At t=0 two electrons are localized at the donor. Two excess 
electrons on AT- AT are self-trapped and not displayed. 



fer occurs activationless for two holes in a GC-GC dimcr 
because U c s — E a \. The estimate of the transfer rate via 
Eq. (|13p at room temperature translates into an ultra-fast 
characteristic time scale of about lOfs. For two electrons 
on the GC-GC dimer, however, the bare U is smaller, 
and the condition U e ff < E a \ results in a smaller rate. 
The transfer rates for two electrons or holes on an AT- AT 
dimer are exponentially suppressed and several orders of 
magnitude smaller than those on a GC-GC dimer. 

At low temperatures the mean-field decoupling of the 
bath modes is not applicable any longer. The quantum 
mechanical nature of bath modes must be taken into ac- 
count, and nuclear tunneling is important: now coherent 
collective oscillations of the electrons become possible. 
As a consequence, the transfer characteristics is changed, 
and the rates can become larger than in the thermal ac- 
tivation regime. 

We employ the TD-NRG to calculate the population 
probabilities do(t) and noA{t) at a temperature T = 
1.5 • 10~ 8 eV for the parameters of the different base pairs 
coupled to an ohmic bath. The real-time dynamics for 
short-time scales is plotted in Fig. [6l The transfer shows 
an incoherent behavior with a few remains of oscillations. 
The decay rates are evaluated by fitting the result to the 
kinetic equations Eq. (fT5|) . These fitted rates are added 
to Fig. [5] as (red) points and are also depicted in Fig. El 

The largest rate is found for GC-GC dimer with two 
electrons, koA — 0.167 • 10 15 l/s, and the equilibrium 
is reached after about 40 fs, while the rate for GC-GC 
(2h) is 25 times smaller corresponding to a characteristic 
time scale of 146fs. For T — > 0, the optimal regime is 
already reached for U c s < E a \ and, therefore, earlier 



than predicted by the Marcus theory. Such behavior has 
also been reported in the biased spin-boson model at low 
temperatures.— Nevertheless, the Marcus rates and the 
TD-NRG rates are of the same order of magnitude in the 
optimal regime. Using the classification of the Marcus 
theory, the transfer in GC-GC (2e) then would be already 
in the inverted region at low temperatures. 

We would like to emphasize that these transfer rates 
are significantly larger than those reported for the single 
charge subspace. The incoherent relaxation occurs on 
time scales between a few fs and 50fs while the envelope 
function of the single charge dynamics decays on times 
scales larger than 200fs as can be seen in Fig. |4] 

For AT- AT (2e) the two electrons on the donor are self- 
trapped in an ohmic bath. For AT- AT (2h) one hole is 
transferred with a characteristic time scale of 146fs. The 
transfer is much faster than in the thermal activation 
regime where the rate is exponentially suppressed. 

A summary of the TD-NRG results is included in 
Fig. El We can identify three regions. Below U c g = 0.5eV 
(left of the dotted line) the system is localized^ for an 
ohmic bath, and the charges are self-trapped. Right of 
the dotted line - white region - we found sequential trans- 
fer. With increasing U e s more coherent transfer is ob- 
served, and the fraction of pair transfer is increasing. 

At low temperatures pair transfer would only be 
present when U c q exceeds AeV for the given reorgani- 
zation energy E a \ in an ohmic bath. The charges remain 
trapped at the donor in the AT-AT (2e) dimer coupled 
to an ohmic bath since U e s = —leV lies in the localized 
regime. This self-trapping is lifted in a super-ohmic bath: 
pair transfer is observed in super-ohmic baths as shown 
in Fig. M 

B. Finite bond-bond interaction 
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FIG. 7: Schematic potential surface E(x) of the lowest states 
for a single bosonic mode as a function of the displacement x 
in arbitrary units for a finite \W\ > t/eff: for W < (upper 
panel) and W > (lower panel). In the latter case *I/t is the 
state with the lowest energy. The hopping matrix element 
between this state and the initial state | t4-, 0) is zero. 



The bond-bond interaction has two effects: First it in- 
duces an effective pair hopping between | fi, 0) and |0, t4) 
and second it lifts the degeneracy between the triplet and 
singlet states |\l/ t ) and \^ s ), respectively. The triplet 
state, which is decoupled from the other states, becomes 
the ground state for a finite W > and U e s > 0. Since 
spin conservation prohibits its population when starting 
initially with two fermions on the donor site, the transfer 
involves only the covalent singlet state Its dynam- 

ics is governed by the energy difference AE — U e s — W 
between |t4-,0) and |* a ). 

For a small U cS = O.leV as in AT-AT(2h) the rate 
predicted by the Marcus theory decreases with increas- 
ing W because the energy difference changes its sign from 
UcH — W > to U e ff — W < 0. These two regimes are 
shown schematically in Fig. [7]for \ W\ > U c s- The depen- 
dency of the rate koA on W obtained at low temperatures 
using the TD-NRG follows the prediction of the Marcus 
theory as shown in Fig. [5] as (red) circles. 

As discussed above, the GC-GC (2h) dimer is located 
in the inverted region as a consequence of the large 



U e g — 2eV. An increasing W reduces the energy dif- 
ference between the initial state and |\P S ) leading to an 
increase of koA as depicted in Fig. [8] (squares). The 
transfer in GC-GC (2e) shows almost no dependence on 
W (Fig. [51 triangles) since it is located in the activation- 
less regime. 

An electron-pair in an AT-AT dimer remains self- 
trapped for \W\ < 0.25eV in an ohmic bath. In a super- 
ohmic bath (s > 1), however, coherent transfer is found 
at low temperatures. The electrons oscillate coherently 
as pairs between the donor and the acceptor while the 
high lying singlet covalent state |\& s ) is only virtually oc- 
cupied. The oscillation frequency is approximately given 

4A 2 

by 77r — fi T777 + 2W c ff where W c s is the renormalizcd pair- 

\u e {{—w\ 

hopping in the presence of the bosonic bath: this fre- 
quency increases with increasing W. 

The occupation probabilities show slow oscillations 
with a period of the order of 1 — 4 ns as depicted in 
Fig.lHl Note that the time scale is multiple orders of mag- 
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FIG. 8: (Color online) Low temperature rate ho a that at least 
one charge is migrating from D to A for AT-AT(2h) (red cir- 
cles), GC-GC (2e) (blue triangles) and GC-GC (2 holes) (yel- 
low squares). The time-dependent population probabilities 
are calculated with the TD-NRG using an ohmic bath and 
a — 0.5. As the transfer is in its incoherent regime the rates 
are obtained by fitting the kinetic equations eq. (fT5|) . 




t/ns 

FIG. 9: (Color online) Low-temperature population proba- 
bility P(t) = do(t) of the doubly occupied donor (starts at 
one) and acceptor P(t) = d,A{t) of the covalent states (starts 
at zero) for W = (blue dashed line), W = O.leV (red dot- 
ted line) and W = — G.leV (black solid line) as a function of 
time for AT-AT (2e) considering a super-ohmic bath (s = 2) 
and a — 1. The population probability of the covalent states 
P(t) = riDA(t) is zero (green squares). 



nitude larger compared to the fs time-scales of the GC- 
GC dimer. The relaxation rate of the envelope function 
corresponds to a significantly longer time scale ~ 50ns. 



V. SUMMARY 

We have calculated the real-time dynamics of a re- 
dox system with multiple charges. At temperatures 
lower than the activation energy we have used the time- 
dependent numerical renormalization group approach; at 
temperatures larger than the activation energy when the 
transfer is governed by thermal activation we have em- 
ployed the Marcus theory. Four different DNA dimers 
have been considered: GC-GC and AT-AT with initially 
two excess electrons or holes on the donor site. All re- 
sults are obtained for realistic matrix elements obtained 
by ab-initio methods.— 

The charge transfer characteristics of a system with 
two excess charges deviates substantially from those of a 
single excess charge. The transfer occurs correlated and 
depends strongly on diagonal and off-diagonal Coulomb 
matrix elements. We predict that an additional charge in 
DNA sensors or DNA devices will show base-pair specific 
signal amplification or attenuation. 

We have calculated the rate k^A that at least one 
charge is transferred from the donor to the acceptor. The 
time-dependent population probabilities of the singlet co- 
valent state (l^s)) and two ionic states (|tl, 0), |0,ti)) 
depend on the difference between the effective diagonal 
Coulomb interaction U e ff and the off-diagonal bond-bond 
interaction W. 

We have found a qualitative agreement between the 
Marcus theory a high-temperatures and the TD-NRG 
at low temperatures. The transfer rate kuA shows a 
dome-shape as function of U e s in the Marcus theory as 
depicted in Fig. [5] It is slightly shifted to lower val- 
ues of U e ff and broadened at low temperatures as calcu- 
lated via TD-NRG. However, self-trapping is only found 
in the TD-NRG. The correct application of the Marcus 
theory requires the restriction to transitions compatible 
with spin-conservation: the triplet state \^t) remains in- 
accessible since we always start with two charges in a 
spin-singlet state on the donor. 

The maximal transfer rate koAt defining the 
activation-less regime, is significantly larger in the sub- 
space of two excess charges compared to the subspace of 
a single excess charge. At low temperatures, the largest 
rate is found for two excess electrons in GC-GC and for 
two excess-holes at room temperature. Electron-pairs 
initially on the donor in an AT-AT dimer are either self- 
trapped (ohmic bath) or are slowly oscillating as pairs 
between D and A (super-ohmic bath) at low tempera- 
tures. 

The dependency of the transfer rate on a finite but 
small bond-bond interaction W can be understood in 
terms of a simple shift of U c s — > U e g — W: The rate 
koA now peaks at s» U c s — W — E a \. Consequently, the 
transfer rate increases for two holes on GC-GC which 
is located in the inverted regime while koA decreases 
for two holes on AT-AT with increasing W. It remains 
constant for two electrons on GC-GC since the rate ap- 
pears to be rather insensitive to W in the activation-less 
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regime. 

Furthermore, the presence of the additional bond-bond 
interaction W splits the degeneracy of the singlet and 
the triplet state, \^> s ) and \^t), respectively. Without 
any further spin relaxation process, the system cannot 
reach the electronic thermodynamic ground state |^t) 
for W > and U e g > since the triplet state remains 
unpopulated at all times. The system remains trapped 
in an intermediate metastable state whose life-time will 
depend on the spin relaxation mechanism not considered 
here. 

The transfer characteristics of two charges in a DNA 
dimer depends strongly on the base pair specific Coulomb 
interactions. These effects will become important when 



such redox systems are weakly contacted. They strongly 
influence the current-voltage characteristics once mul- 
tiple excess charges are relevant at higher source-drain 
voltages^ 
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